CONSTRUCTION OF THE THERMODYNAMIC JAMMING 
LIMIT FOR THE PARKING PROCESS AND OTHER 
EXCLUSION SCHEMES ON Z'^ 



THOMAS LOGAN RITCHIE 

Abstract. We provide an explicit construction of the thermodynamic jam- 
ming Hmit for the parking process and other finite range exclusion schemes on 
TJ^ . By means thereof, a strong law of large numbers for occupation densities 
is accomplished, and, amongst other results, the so called "super-exponential" 
(i.e. gamma) decay of pair correlation functions is established. 



1. Introduction 

In its simplest form the car-parking process on Z** (also called random sequential 
adsorption (RSA) in physical sciences) may be (informally) described as follows: 

a site is chosen uniformly over the box A^f'' — {~n, . . . ,0, . . . ,71}'^ (the target re- 
gion) and is occupied (spin/state 1) provided all its nearest neighbours are vacant 
(spin/state 0). Once occupied a site remains so for ever. The process continues 
until all sites in h!:h^ are either occupied or have at least one of their nearest neigh- 
bours occupied (the so called jamming limitj. Call Pn the terminal configuration of 

This model as many of its generalizations (including longer range interactions and 
particles of different kinds) have long received attention in physical sciences, wherein 
issues like kinetics near the jamming limit, spatial correlations, assymptotics (as the 
target region increases to infinity) of the first moments of |P„| (the ultimate num- 
ber of adsopted particles, i.e. occupied sites), amongst others, have been studied 
by a series of non-rigorous theoretical methods, computer simulations and experi- 
ments piTT]. 

Until recently, mathematical results have been restricted to one-dimensional mod- 
els [3 [121 HHIini [1113 and to quasi-one-dimensional models [7|[T1[2]. 

In 2002 M. D. Penrose ^B] proved, amongst other results, a law of large numbers 

(LP -convergence) for the sequence ( |-Pn|/|An''''| J for a large class of exclusion 
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models (including those with finitely many spin states) , treated therein as (coupled) 
finite interacting particle systems. 

In this paper we construct in a common probability space (i.e. couple) a limit 
"parking" algorithm P, having the whole Z"^ lattice as target region, together with 
the finite parking algorithms P„, n > in such a way that P„ — > P almost surely. 
This construction being sufficiently general to cope with all sort of finite range 
(exclusion) interactions in a {0, 1} spin state scenary. In this set up P may be in- 
terpreted as the jamming limit of a virtual parking process over the whole 'Z'^ lattice. 
We call P the thermodynamic limit of the sequence {Pn)n>Oj and its distribution 
fj, the corresponding thermodynamic limit measure. 

By means of this construction (section |3Jl , we prove a strong law of large num- 
bers for the sequence ( |P„|/|a1''-' I J , via the classical multiparameter ergodic 

V / n>0 

theorem (section 0J; the super-exponential (i.e. gamma) decay of pair correlation 
functions fsubsection l5.1(l : the super-exponential speed of convergence^ to the ther- 
modynamic limit measure /i fsubsection l5.2|l and discuss how to obtain lower and 

upper estimates for p, the almost sure limit of ( |P„|/|Ai^''| ) (sectional. 

V / n>0 



2. Basic Definitions and Notation 

The following structures will be basic blocks in the forthcoming constructions. 
Therefore we take the opportunity to establish the notation to be used henceforth. 

Probability and Topological Spaces. The random objects we shall deal with 
in this paper will be defined on a fundamental probability space (i) and assume 
values in a Polish measure space^ (ii). As construction is the realm of this paper, 
we make these spaces explicit: 

(i) : (^{0, 1)^'' , X^Q ^-^^ , where S(o,i) is the usual Borel cr-algebra in (0, 1), 

A(o,i) is the Lebesgue (uniform) measure in (0, 1) and means product in 
the appropiate sense; 

(ii) : ^{0, 1}^'', yB^Q , where ^{o.i} is taken for the power set of {0, 1}, 

7'({O,1})'^='{{O,1},{O},{1},0}. 

Since both (0, 1) and {0, 1} are Polish topological spaces (under the euclidean 
and discrete metric respectively), so are (0, 1)^ and {0, 1}^ (under the respective 
product topologies). Moreover the product cr-algebras S^g^^ and are the 

respective Borel cr-algebras of these Polish spaces. Therefore the measure spaces 
above, viz ^(0, 1)^'' , B^^ ^-^^ and ^{0, l}^'' ,B^q , are naturally related to the Pol- 
ish spaces (0, 1)^ and {0, 1}^ - product topologies implicit- in the sense that the 
(T-algebras of the former are the Borel cr-algebras of the latter. 

These remarks may seem pedantic at first glance, but they are the key to conclude 
convergence of distribution from (almost sure) pointwise convergence of general 
random objects, as we shall see below. 



^In a sense to be specified in subsection 15.21 

■^a measure space whose cr-algebra is generated by a topology which is metrizable, complete 
and separable. 
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Graphs. In what follows the notion of vicinity between sites of the lattice will 
be implemented by the graph =^ ^Z**, E^^^ , where Si'''^ is taken for set of bonds 

joining distinct sites x,y in Z^, such that \x — y|sup '= maxi<j<(i \xi — yi\ < v, for 
some distance z/ e N* N \ {0}. 

Parking Schemes. A subset si"^^ C {0, l}'^-"' , where A^''^ {-v, . . . ,0, . . . ,1^^ 
is the box/ball with centre at the origin and radius u in the Z** lattice, such that 
si^^ 3 O"^^**' will be called a parking scheme (with radius u of interaction) on the 

Z'^ lattice.'' 

Parking Processes. For a given parking schem,e s'if \ we define the parking pro- 
cess on the box A^'^^ with null boundary condition, PPnbc, to be the 
following (random) algorithm: 
Step 1: set P„ = QF"; 

Step 2: choose a site x G uniformly amongst the sites not chosen previ- 
ously; 

Step 3: if 0^ (P„)| A^,'^) e sif \ then P„(x) ^ 1; 

Step 4: if there are points in not chosen yet, then go back to step 2, else 
stop the algorithm. 

Above, for all y € Z'^, 9y is the translation operator, i.e. for all P G {0, 1}^ , 
0y (P) (x) '= P{x + y). And P| aS'*'' means the restriction of P to the box A^,'''*, 

that is, P| A^"^) e {0, 1}^^"' and for all x G A*.'^^ P| Ai^\x) = P{x). 

It should be clear that the algorithm described above generates a random element 
P„ in the PoHsh measure space ^{0, l}^**, B^g and therefore a corresponding 
probability measure is induced therein. 

We put these ideas into a probability setting by means of 

Definition 2.1. For a given parking scheme s'if\ we define the parking process on 
the box An^ with null boundary condition, PPnbc, to be the following 

algorithm: 

For eachu e {0,lf\ 

Step 1: set P„(cj) = 0^" ; 

Step 2: choose x G A^'''' such that 

u>{x) = inf{a;(^) : ^ G A^f^ and ^ has not been chosen previously}; 
Step 3: if 0^ (P„(a;))| A^''^ G si'^\ then P„(a;)(x) ^ 1; 
Step 4: if there are points in A^^ not chosen yet, then go back to step 2, else 

stop the algorithm. 

□ 

The formal algorithm described above is probabilisticaly equivalent to the former 
(informal) algorithm in the sense that both induce the same probability measure 



^As it will become clear from the forthcoming definitions, we may suppose, without loss of 
generality, that r]{0'^) = for all r] e sifK 

Henceforth we shall abuse notation and write for 0**, whenever no confusion may arise. 
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Hn in the Polish measure space (^{0, 1}^ jl^fo i})- This equivalence being due to 
the fact that the minimum of an iid finite family of uniform random variables is 
uniformly distributed amongst them. 

We call the output Pnito) the jamming limit of the parking process just de- 
scribed^ and observe that P„ is a {0, l}^''-valued random element on ^(0, 1)'^'^ , B^^ -^^^ 

' -^(0,1)^- 

Finally, we remark that each Pmn > 0, is well (i.e. uniquely) defined, except for 
a set of null probability; ambiguity might occur at cj's such that lu{x) — uj^y) for 
at least two distinct sites x,y S Z'^. We shall return to this point later on, but, as 
it usually happens in probability theory, that should pose us no difficulty at all. 

A natural question that arises at this point is the existence (or not) of a ther- 
modynamic limit probability measure /i for the sequence (/i„)n>i of probability 

measures on (^{0,1}^'' ,B^^ induced by the sequence {Pn)n>i of random ele- 
ments just defined. We shall delve into this question in the next section. 

3. Construction of the Thermodynamic Limit 

In what follows we show that, for a given parking scheme sif'\ there exists 
a unique probabihty measure on ^{0, 1}^'', such that ^„ ^. ^ We 
shall accomplish this task by means of constructing a random element/algorithm 
P : ((0, l)^^ef;l), Agi)) ^ ({0, lV\Bf; ,^) such that P„ -> P a.s. ' 

For that purpose we set the following definitions: 

Definition 3.1. A path of points in the graph Q^f"^ is a finite sequence (2:j)o<i<n' 
Xi e Z'*; such that \xi+i — a:^j|sup <i^,0<i<n — 1. □ 

In what follows we shall say that the path (a;i)o<i<n has length n (rather than 
n + 1), bearing in mind that n bonds participate in such a path. 

Definition 3.2. A path (a:^i)o<i<n ^'^^'^ decreasing (subject to ^^-"£(0, 1)^ ), 

whenever the numerical sequence (iCi(w))g<j<„ is (strictly) decreasing.^ In this case 
we write {^i)o<i<n' "'^ simply | (a;i)o<j<„, should u) be clear from the context. 
□ 

It should be clear from definition 13.21 that a decreasing path is always self- 
avoiding. 

Definition 3.3. For x,y & and lo G (0, 1)^ , we say that x is influenced by y 
(or, equivalentely, that y influences x) subject to uj, if there exists a decreasing path 
{^i)o<i<n' such that Xq — x and Xn = y. In this case we write y t(t^) x, or 
simply y 'I X, should co be clear from the context. □ 

We remark that the f relation is both transitive and reflexive. 



^In the forthcoming text we shall frequently abuse torminilogy and use P„ to denote also the 
underlying parking process leading to it. 

^the symbol "=>" denotes convergence in distribution. 

'^pointwise convergence regarding the product topology in {0, 1}^ . 

^Xi(Lo) is taken for Lj(xi), which should be formally prefered. 
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Definition 3.4. Given a subset X C 'Z,'^ , we define the armour of X (subject to 
LO G (0, 1)^ ) to be the (random) subset 



□ 



As before, we shall omit lo whenever no confusion may arise. 
It is clear from definition 1231 that, for X,Y dZ'^ and uj e (0, 

(i) : A{X){lu)dX; 

(ii) : X = 1±J:Li^.^^WH = U':Li-4(X,)H; ' 

(iii) : X CY^ A{X){uj) C A{Y){uj). 
What may not be so clear is the following lemma. 

Lemma 3.1. For any finite X C Z'^, A{X) is (almost surely) finite. 

Proof: 

We divide the proof into three parts: 
Part (i) 

Suppose initially that X = {0} and consider the event An (^A{X) (t aI'^^^, 

n G N. An amounts to saying that there exists y G ifi \ hik} such that y t Oi 
that is, there exists a decreasing (and hence self-avoiding) path (xi)o<i<m, m > n, 
starting at and ending outside the box A^fJ . Thus An occurs only if there exists 
a decreasing path of length n, {xi)o<i<n starting at 0. Now observe that the 
probability that an arbitrary (but fixed) self-avoiding path (xi)o<i<n is decreasing 
is l/(n -I- 1)! and that the total number of self-avoinding paths starting at with 
length n is no larger than (2i^ + l)"* ", to conclude (by subaditivity) that 

, , (2j^+ I)'''" (2iy4- 1)'*-" 



(n-l-1)! ~ n\ 

Since the sequence (A„)„>o is decreasing and P(^,i) — > by (^, we must have 
P(n^=o^") = 0- Now ob~serve that {^A{X) = 00} = CC^o^n to conclude that 
P({#yt(A) < 00} = {#A{X) = 00}") = 1, that is A{{0}) is almost surely finite. 

Part (ii) 

Set X = {x}, for some x G Z'^. 

In this case we observe that {#yl(X) < 00} — 9^^ {{#A ({0}) < 00}) and make use 
of the stationarity of the product measure A^p to conclude that P {{4j^A{X) < 00}) 
— P ({#-^ ({0}) < 00}) = 1? by part (i) above, that is A {{x}) is almost surely finite 
for aU xeZ"^. 

Remark: we might have observed the event equality 

Ai{x}) <t B(x,nv) = (^({0}) AW) 
made use of stationarity to yield 

P(^(M) <t B{x,niy)) = P (^({0}) <t Ai^)) = P (^({0}) A(5) 



'l+J denotes disjoint union. 
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and followed the same steps of part (i) to see that A {{x}) is almost surely finite. 
This translation technique will be recurrent in the forthcoming proofs. 

Part (iii) 

Suppose X is a finite set, that is X — l+J"^]^{a;i : Xi E X}. Thus by the remark after 
definition E3|^(X) = [Jti^ii^^}) ^ {M^) = 00} = ULi {^^{{x,}) = 00} 
^V{{#AiX) = 00}) < Er=i IP ({#-4 ({a: J) = 00}) = by part (ii) above. That 
is A{X) is almost surely finite. 

■ 

In what follows it will be convinient to generalize parking processes on boxes (defini- 
tion!^^ to parking processes on arbitrary finite subsets of Z''. This generalization 
does not involve any difficulty and is performed in the natural way. 

Definition 3.5. Let X be a finite subset ofU^. For a given parking scheme si'^\ 
we define the parking process on X with null boundary condition, 

(^X Sjf"^^ PPnbc, to be the following algorithm: 

For eachuje (0, 

Step 1: set Px{oj) = 0^^■ 
Step 2: choose x E X such that 

Lu{x) = inf{ix'(^) : ^ G Xand ^ has not been chosen previously} ; 
Step 3: if 9^ {Px{oj))\ t^^f' € s'f^ , then Px{'^)ix) ^ \; 

Step 4: if there are points in X not chosen yet, then go back to step 2, else 
stop the algorithm. 

□ 

As before, we call the output Px (w) the jamming limit of the parking process just 
defined and observe that Px is a {0, l}^'*-valued random element on ^(0, 1)^'' ,B^q 



Now we are ready for the definition of the limit random element/algorithm P: 
Definition 3.6. The limit random element/algorithm 
P: ((0,l)^^6g,),Ag,)) - [{0,1}^\B%,^) ^s defined by 

Picj)ix) F4({.})(.)(u;)(x), Vx e 

□ 

In words definition 13. (jl above means that for each u> e (0,1)^'' and x € Z*^, 
we determine the armour of the set {x), A{{x]'), by means of definition 13.41 
construct/implement the parking process over ^({cc}) according to definition 13.51 
and then assign the value of the corresponding jamming limit Pa{{x)) site x to 
P(uj)(x). Of course a parking scheme S^f' is implicit in definition 13. 61 

Our next step is to prove that Pn P a.s. For that purpose we shall need the 
following lemma: 

Lemma 3.2. Pn\A{{x}) = Pa{{x}}\ A {{x}) in {A„ D A{{x})}. 

In words lemma rOl asserts that the (random) sets A {{x}) are stable with respect 
to the sequence {Pn)n>v sense that, once n becomes sufficiently large so that 

A„ D A {{x}), the jamming limits Pm, m> n are identical over A {{x}). 
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In the proof of lemma 13.21 the notions of inner and outer boundaries of a set 
X (lU^ will be employed. We make them clear by means of the next definition. 

Definition 3.7. For each i' >1, the ly— inner boundary of X CZ 'Z'^ is defined to be 
the set 

(5f,™'(X) =^ {x eZ'^ : X G X and 3y e Z'^ \ X such that \x - y|sup < i^} 
In the same fashion, the v-outer boundary of X C Z'^ is defined to be the set 
{x eZ'^ : X i X and 3y e X such that \x - y\snp < i'} 

□ 

For simplicity the subscrits will be suppressed, whenever v is understood from 
the context. 

Remark: an immediate consequence of definitions 13.41 and 13.71 is that a;(a;) < 
y{uj), whenever x € {A{X)), y € (5(°"*) {A{X)) and \x - y\sup < v, for all 
X C Z''. That is, an armour has the noticeable property that in the interface with 
the outer enviroment its sites are numerically smaller than the interacting outer 
sites. 

Proof of lemma 13.21 Initially we observe that, for x e Z**, (almost all) u G 
(0, 1)^ induces a natural order in the set ^({x}) (w) C Z'', in the sense that, 
for y,z £ A {{x}), y ^ z, whenever y{uj) < z{uj). 

Thus, for each uj, one can construct a sequence iyi)i<i<#A{{x}){uj) in -^({2;}) 
such that yj =4 Vk, whenever j < k; this sequence being (almost surely) uniquely 
determined. 

Now for a fixed u! and a time r G {1,2, . . . ,^A{{x}) (ui)}, we observe the 
parking processes Pn (definition and Pa{{x})(uj) ^definition I3.5|l at time r, i.e. 
immediately before an arrival (step 2 of the corresponding algorithms) at site yr- 

We conclude the proof by induction (over t): 
Base of Induction: 

At time r = 1, Pn\A{{x}) = Pa{{x})\A {{x}) = Q-^d^}) , by step 1 of respective 
definitions. 

Inductive Hypothesis: 

Suppose that Pn\A{{x}) = PAi{x})\A {{x}) at time r G {1,2, . . .,4A{{x})}. 
Inductive Step: 

If yr 4- i5i>"'^ (^({3^})) 5 by the inductive hypothesis, Pn and Pa({x}) are identical at 
time T in the (^-neighbourhood of yr, i.e. B{yT,i') {y £ Z'^ : \y — t/rlsup < i^}, 
hence 0j,^(P„)|A[.''^ = Oy^ 

{PA{{x}))\^^y''' ■ Thus, by step 3 of the corresponding 
definitions, we see that F„|yl({a;}) = ^U({z:})|.A ({cc}) at time r + 1. 

If yr € i^t'*"'' (^({a^}))i again by the inductive hypothesis, Pn and Pa({x}) coincide in 

B{yr,v)C\A{{x}). But, since B{yr,v)C^{A{{x})Y C 5i°'"^\A{{x})) , no P„-arrival 
may have occured in this region up to time r by the remark after definition 13.71 
and therefore both P„ and Pa({x}) must equal thereon. Hence P„ and Pa{{x}) 
coincide again on the whole B{yT, v) at time r, and by the same reason as before 
we must have P„|^({a;}) = Pa{{x})\A{{x}) at time r + 1. ■ 
Lemma |3. 21 is the key to the main result of this section: 
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Theorem 3.3 (Almost sure convergence to the limit algorithm P). 

Pn P ^fo,l)'^-^- 

Remark: before proving theorem 13.31 we emphasize that the a.s. convergence 
therein is relative to the usual product topology of {0, 1}^ . 

Proof: 

Bearing in mind the foregoing remark, it is enough to show that for A^q -^-j -almost 

all uj G (0, 1)^ and every to > 0, P„([x))|Am'' = P{'-^)\-^rn , for n sufficiently large. 

Now, given to > 0, take n large enough so that AiAin') C Aif\ This is possible 
for almost all u> in view of lemma l3.1l For each x £ Am\ A {{x}) C A{Am^) (remark 
after definition 13. 4|l . Hence 

Pix) '^'^P^a^y^ix) '^"J^p„(x), vx e A(;J) 

and we are done. I 

Since the random elements (P„)n>o and P take values in the Polish measure space 
^{0, 1}^ ,K|Q j^j^ , almost sure convergence implies weak convergence of the corre- 
sponding distributions therein (111 prop. 9 pag.355]. Thus the following corollary: 

Corollary 3.4. // (/in)n>o ^.i^-d fJ. are the corresponding probability distributions of 
{Pn)n>o and P in (^{0, 1}^'' , B^^ , then ^„ =4> fi. 

Before finishing this section we wish to emphasize that the construction of the 
limit measure carried out in definition 13.61 endows it with two noticeable properties, 
viz stationarity (i.e. translation invariace) and ergodicity. 

In order to formalize this point a little further we observe that the random 
elements (P„)„>o and P just defined f definitions 12 . II and 13 . 51 respectivelv) are well 
defined, except for a set of probability zero, viz the set of w's in (0, 1)^ that assign 
a same numerical value to at least two distinct sites of Z"* or on which an infinite 
descending path occurs. 

This difficulty may be easily removed by setting P„ = P = 0^ , n > 0, for 
instance, over this troublesome set. Once this is done, we can define the real 

measurable function / on ^(0, 1)^'' ,B'^^ ^fo*!)) /(^) Pi^)i^) ^^'^ observe 

that P{uj){x) — f {9x{uj)), for all x in Z"^. Since the product measure A^q is 

ergodic and stationary in ^(0, l)^'' , B^^ , so will be /i, the distribution of P 



m 



({o,iF 



'^{0 i}j ■ This is a standard argument in ergodic theory and we refer the 
reader to |31 chapt.6] for further details. 

Therefore we have established another corollary to theorem l3.3l 

Corollary 3.5. The thermodynamic limit measure /i is both stationary and ergodic 

The construction of the thermodynamic limit measure fj, carried out in this sec- 
tion corresponds to the concept of perfect simulation in the sense of Fernandez, 
Ferrari and Garcia In this technique, with the aid of a finite (yet random) 



number of i.i.d. (uniform) random variables, one generates (i.e. simulates) accord- 
ing to a prescribed thermodynamic measure a finite portion (window) of a infinite 
random object (field). 

Besides the straightforward advantages comparing to ordinary Monte Carlo sim- 
ulations, wherein only finite approximations of the thermodynamic measure are ac- 
tually simulated, this constructive technique is also very useful regarding theoretical 
proofs as the forthcoming sections testify. 

4. A Strong Law for Occupation Densities 
First of all we formalize the concept of occupation density refered to in section^ 

Definition 4.1. For n G N, we define the occupation density at the jamming limit 

Pn to he the random variable 

def S.gAW Pn{^){x) / , , , 

□ 

As we mentioned in section^ the existence (or not) of lim„^oo lE(p„) is a ques- 
tion of relevance in random sequential adsorption (RSA) contexts. Untill recently an 
(affirmative) answer was known only for one-dimensional models [Ul I12L [T^ I14L [T^ 
or for quasi-one-dimensional models Ql [3 ■ Recently Penrose ^3 proved a law 
of large numbers (convergence in probability) for the sequence (pn)ri>i for a large 
class of exclusion schemes (considered therein as finite particle systems) on Z'^. 
This class overlaps to some extend to the class considered here.^° Naturally, this 
result assigns an affirmative answer (by means of a straightforward application of 
the Bounded Convergence Theorem) to the existence of lim„^oo E(p„). 

Here with the aid of the construction carried out in section|31 we shall accomplish 
a simple proof for a strong law of large numbers (SLLN) for (p„)„>i. 

The first step to this proof is to observe that we can invoke Wiener's multipa- 
rameter ergodic theorem jlO. Appendix 14. A] to conclude that 

lim ^"g^"' ^ ' = lim ^-^f-'' p ' ^ E(/) 

n^oo \K^n \ (2^+1)'' 

(2) = E(/o0o)=E(P(O)) Afo'i)-a.s. 

for some constant p e [0,1]. This is true since / is measurable and integrable 

on ^(0,1)^ ''^fo 1)' '^fo 1)) ^'^^ ^^'^ product measure is stationary and ergodic in 

((o,lf^<,r))• 

One might guess that the discrepancies between P{uj)\K^'' and Pn(u})\h!'n^ are 
due to some kind of boundary effect. Next lemma quantifies the extension of this 
boundary pertubation. 

Lemma 4.1. If An is the event ^aI'*'^ ^ a|^^^^^j^^|, then 

P(limsupA„) — 0, i.e. with probability one the sequence (A„)„>i occurs finitely 
many times.^^ 

-'^'^We refer the reader to section [5] for a discussion about this point, 
greatest integer less than or equal to x. 



Proof: 

By means of the estimation carried out in lemma IXTI and bearing in mind the 
stationarity of the product measure as we emphasized in the remark after part (ii) 
of lemma ISm subaditivity yields: 



P(A„) < (271+1) 

Therefore by Stirling's formula^^: 



bn 



But, since 6„ ~ o(l/n ), 



2\ 13 



oo 

On < OO 
n— 1 n— 1 



< OO. 



6„ < OO 

n=l 

Thus the proof is concluded by an application of the first lemma of Borel-Cantelli. 



Now we are ready for the main result of this section: 
Theorem 4.2. A Strong Law of Large Numbers for Occupation Densities: 



Pn 



P 



(0,1) 



-a.s. 



Proof: 



We know by lemma UTI that the last occurrence of {An)n>i is almost surely finite, 
that is, n = sup{n e N : A„ occurs} < oo a.s. Given e > 0, pick i e N sufficiently 
large so that i £ [n + [-^nji^, (n + 1) + + for some n > n and such that 

|^ 2(r»+i)+|vm i:2^^'^_ ^ ^ > Therefore PlAi"^) ^ P.IAI'^' by the 

proof of theorem 1^21 
Hence 



{2i + l)'i 



< 



{[(ti+1)+LV^J->]-2+i}''-(2«+1)'^ 
(2ri+l)'' 



That is 



Since lim, 



2(K+l)+[Vn+lJ -21^+1 
(2n+l) 



hm — . ' , , J T— — -— 1 — = a.s. 



1 < e 



(2i + 1)'^ 



(2i + 1)'^ 



(24l)'^ ^ P ^y®, P. = TmTTvJ ^ P a-s- 



(21+1)" 



As we mentioned in the beginning of this section, a straightfoward application of 
the Bounded Convergence Theorem now garantees the existence of lim„^oo E(p„). 

^^n! ~ n"e~"y'2T7n, i.e. ri"e~"\/27m/n! — > 1 as n — > oo. 
"'^'^Landau's notation: /„ = o(g,i) means that f,i/gn — > as n — » oo. 
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Corollary 4.3. lim„^oo IE(p„) = p 

In section El we briefly discuss how to estimate the constant p numerically. 

5. Super-exponential Behaviour 

First of all we precise the meaning of super- exponential decay. 

Definition 5.1. A function / ; N — > M is said to decay super-exponentially when- 
ever f{n) = o(e-""),Va > 0, i.e. lim„^oo /H-e"" = 0,Va > 0. □ 

In physicochemical literature the so called random sequential adsorptian (RSA) 
processes are distinguished from equilibrium systems by several noticeable proper- 
ties |H1 El ■ Amongst them we pay special attention to super-exponential decay of 
pair correlation functions, a feature of RSA models that differs remarkably from the 
(usual) exponential decay of equilibrium systems. Although described in physical 
and physicochemical literature |^ I17|. a rigorous account thereof is available only 
for one-dimensional models T. 

Here super-exponential decay of pair-correlation functions is established from 
the construction of the thermodynamical limit measure p, carried out in sectional 
This amounts to saying that this feature is a universal property of the whole class 
of exclusion schemes studied in this paper, regardless of dimension (c?) and range 
of interaction (j/). 

Another universal property derived from section|31is the super-exponential speed 
of convergence to the limit measure p. We cope with this point in subsection 15. 21 

5.1. Super-Exponential Decay of Pair-correlation Functions. We recall that 
the correlation coefficient <tx,y of two random variables AT, Y is defined by ax,Y = 
Cov{X, Y)/ (axo-y), where Cov{X,Y) stands for the covariance of X and Y, and 
CTx, o'Y for the corresponding standard deviations. 

Regarding the exclusion schemes studied here, one might guess how the state 
of site X G influences the state of another site y € 1/^ with respect to the 
limit measure p. Therefore we shall investigate the assymptotic behaviour of the 
correlation coefficient a,j.y of the random variables P{x) and Py P{y) as 

the distance \x — y\sup tends to infinity. Thanks to the translation invariace of p 
fcorollarv l3.5|l a,j.y = <7o(^y-x), and we may suppose x = without loss of generality. 

The central result of this subsection (theoremjOJ asserts that <tqj. decays super- 
exponentially as x tends towards infinity. More precisely, this amounts to saying 
that limj^^oo croa;e"l^l = 0,Vq; > 0. In order to prove theorem 15.11 it will be con- 
venient to enlarge (i.e. square) our basic probability space (j^O,!)^'' ,8^^ X^^ ^ 

to 1 [(0, 1)^] , j-)2, ) • In this new setting each u) e [(0, 1)^] may be 

thought of as an ordered pair Cu — (wi, 0^2), for uii S (0, 1)^ , i = 1,2. 

On this squared probability space we define the following projection operators 

taking values in the former unsquared measure space (jO, 1)^ ,8^^ 

Definition 5.2. 

Cardinal Projection Operators: 

/ ~ ^ dcf / ~ ^ dcf 

UJl(UJ) — OJi; U!2(0J) = OJ2. 
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Hybrid Projection Operators: 

for n e N and bearing in mind the canonical representation x — {xi,X2, 
. . . , xj), for sites in as d-tuples of integers}'^ □ 

In words definition 15.21 denotes that a;|"'' plays the role of uji in the half hyper- 
space {xi < n} and the role of ujj in {xi > n}. 

A straightfoward implication of definition 15. 21 is that all the projection operators 
above share the same probability distribution in B'^q viz A^q -^y Moreover, wi _L 

UJ2 and Jj^^ ± w^"^,Vn > 0. 

The foregoing paragraph amounts to saying that -P(ij-'i), P(y>2), P{'jJi^I), P{lo^i) 
share the same probability distribution /i in Bf^ -^^y and that ^(^["3) and F(w2"i) 
are independent, i.e. 

(3) ^ P{cj2) ^ P{J"i) ^ P{4"I) ^ M,Vn > 

(4) P{J;;i)±P{4"^),yn>0 

Now we are prepared to prove the key result of this subsection, viz^^ 

Theorem 5.1 (Super-exponential decay of correlation functions). 

a-Qx decays super- exponentially as x —f 00, i.e. lim^;^ e"l^l = 0,Va > 0. 

Proof: 

Let x = {xi,X2, . . . ,Xd) € Z'' and suppose initially that xi — \x\ = supj^^j^^; 
{|xi|}. Set y — [2^1/2] and pick the largest n g N such that 2vn -\- < xi, i.e. 
n = Yxi/2v\ - 2. Hence^'^ nv + 2v <y <\x\- {nv -{- 2v). 

Now observe that A{{Q}){uji) C B{{),nv) => A{{{)}){uji) = A{{Q}){uj'fl) . This is 
true because wi and uj'f^ are indistinguishable on the left of y f definition 15. 2|l . By 
the same token A{{x}){uji) C B{x,nv) => A{{x}){uji) = A{{x}){lu!^^1). Thus on 
the event E e B^^ ^-^2, 

E =^ (^({0})(a;i) C B(0, nv)) n {A{{x}){uoi) C B{x, nv)) 
we must have 

P(a;i)(0) '° =^P^({o})(.oK)(0) = ^^({o})(.S)(-^!'2)(0) '^'-^^'(^S)(0) 

P{u.,){x) = ^^({.»(.<y))('^?i)(^) ^'''^P{^l'l){x) 

and hence 

(5) P{oJi){0) ■ P(^i)(x) = F(^g)(0) • P(a;g)(x) on E. 

^^l{a;i<n} = 1 if 2^1 < and otherwise . . . 

stands for independence. 
^^Henceforth, for x £ JA, \x\ stands for |a:|sup. 

"'^^Remembor that v stands for the radius of interaction of the underlying parking scheme s[f^ . 
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Since 

(6) P{E'=) < 2 

by (UJ; (01 and © yield 



{2iy + I)''-" 



(7) 



E 



< 



However 



(8) 



(|([F(c.i)(0) • P{cu,){x) - pU^^W) ■ P{4'})(x) 



E[p{J^lm ■ Pi4'}){x)] = E[PiJ^})m-E[p{4y}){x)] 

= E[PK)(0)] •E[P(c^i)(a;)] 



where we have used in the former equahty and 121 in the latter. Thus substituting 
© in O yield 

(9) |E[PK)(0) • P{uJi){x)] - E[PK)(0)] • E[P(c.i)(a;)]| < 2 • ^ ^, ' 

Meanwhile, P o uji (on the enlarged/squared probability space) and P (on the 
former/unsquared probability space) share the same distribution /i in B^^ ^ thanks 
to the remark after definition 15. 21 i.e. (Pouii) ~ P ~ /i. Therefore O amounts to 

(2z/+ l)"*-" 



(10) 



|E(Po-P,)-E(Po)-E(P,)| <2 



Substituting n in terms of as in the first paragraph of the current proof, and 
recalling the definition of covariance of two random variables, we finally have 

{2iy+ l)''-(LI^I/2i'J-2) 



(11) 



|Cov(Po,P,)| < 2. 



(L|x|/2H-2)! 

At this point we observe that nothing would change in the foregoing calculations if 
it were the case that Xi = ±\x\, 1 < i < d, in the beginning of this proof. A suitable 
redefinition of the projection operators uj[ \ and uj^ \ along the corresponding coor- 
dinate axes being all we would have to do to establish once again. Therefore 
remains true for all x ^lA. 
Having said this, since 

al =f Var(Po) = Var(P,) =^ a^, V:r G 1'^ 

thanks to the stationarity of /i fcorollarv l3.5|l . we simply make use of the definition 
of correlation coefficient in to establish that 



(12) 



2_ {2v+lY-iVM/^-\-^) 



(e-"l^h, Va > 



and hence conclude the proof. 
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5.2. Speed of Convergence to the Thermodynamic Limit. Since P„ actually 
takes values in the finite set {0, 1}"^" \ it is reasonable to quantify the rate at which 
the sequence {pn)n>o converges to the limit measure /i by estimating how fast the 
numerical sequence IfJ-iL) — /x„(L)| converges to zero, whenever L is a local set, i.e. 
a set determined by a finite number of sites in Z'^. 
The concept of locallity can be formalized as follows: 

Definition 5.3. A set L G j^pd is said to be local (or cylindrical) if for some 
m £ N there exists a subset L C {0, 1}^™' such that 

{0, 1}^' 9 c e L ^ {0, 1}^" 3 oj\Ali^ e L 

We write C for the class of all local sets in B^q . □ 

Next theorem asserts that the numerical convergence above is super- exponen- 
tially fast. 

Theorem 5.2. If L e C, then - M„(i)| = o(e"""), Va > 0. 

Proof: 

Choose L and m like in definition 15. 31 For n > m, 



P [a{A(£) A„) < {2m + l)'^ . ^'^''t}^'!: ^ 
by and subaditivity.^* 



However P\h!£> = P„|A^^^ on {A{h!w,^) C A„j by the proof of theorem lO 

^l{L) = P({P e L}) = P({P|A(,? e L}) = E (i^piaL-)^^}^ 

= P({F„ e L}) = P({P„|AW e L}) = E (l{p„|A(f ^l}' 



Moreover 



and 



Hence 



|/i(L)-/in(i)| 



E 



< 
< 



'(^(A(f)^A„) < (2r.+ l)^.i^^i±g^ 



o(e^""), Va > 



[x] = the least integer greater than or equal to x. 
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6. Rigorous Numerical Estimates and Final Remarks 
We know from section^ that 

lim p,=p = E(Po) = F ({Po = 1}) = 1 - P ({Po - 0}) 

i — >oo 

Moreover, the construction carried out in sectionOhas shown that Pq is determined 
by the relative order of a super-exponentially smaU^^ random number of the i.i.d. 
uniform random variables uJx,x G 7/. Therefore, bearing definitions 13.51 and 13 . 61 in 
mind, it is possible to identify by inspection disjoint subsets of both {Pq = 1} and 
{Po = 0} and hence determine lower and upper estimates for p. 
We clarify these ideas by means of a concrete example: 

Set d = 1, = 1 and si'^^ = S'f^ = {(0,0,0,)} as in section El This set-up 
amounts to the one-dimensional version of the nearest neighbour exclusion scheme 
we touched on in section ^ Define r to be the smallest non-negative integer where 
a local minimum with respect to occur, that is r{uj) '= mi{i G 1^+ : uji < 

Wi+i}. Similarly set l{uj) — sup{i G Z_ : Wi-i > uji}. 

Now consider the partition {Eij; i, j G Z+} of defined by^° 

{uj G (0, 1)^ : l{uj) = i and r{uj) = j} 

and set pij = F{Eij). 
It is easy to see that 

oo 

{Po = 1} = {CJ G (0, If : P(cc>)(G) = 1} = 1+) E2^,2J 

i,j=0 



Hence 



p = P({Po = 1}) = E IP(^2.,2j) - P2^2, 
iJ—0 



Since PiE.j) 



(i + j-K3) • + ■ U + (j + .7+3) ■ (i + j + 2) ■ + - il-j! 

+ (^ + J+3)•(^ + J + 2) ' + + {j +\)\ ■ i\) Vz,jGZ+ 

by elementary combinatorics. It is possible to generate lower bounds for p by 
summing up the first even- even probabilities and to generate upper bounds by 
summing up the first non-even-even probabilities. Thus 

2 

0,4304- P2^ao < P 

hj=0 

/ 



< 1 - 



P2i,2j+1+ Y P2i+l,2j + Y P2i+l,2j + l = 0,4339 



0<i<2 0<i<l 0<j.1<l 

0<j<l 0<j<2 - 



19: 



in the sense of inequality flTl. 
20z+ Z+ U {+00}. 
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It is interesting to observe that these estimates are in agreement with the assymp- 
totic first moment of the number of adsorpted particles calculated by J. K. Macken- 
zie by means of analytical tools such as Laplace transforms, differential equa- 
tions, etc...^^ 

Of course the same readiness is not present in higher dimensional settings. Yet 
it is possible to obtain non trivial estimates for p at least for simple models such 
as the two-dimensional nearest neighbour exclusion scheme and, to our knowledge, 
even in this case no rigorous estimates for p exist. 

Parking Processes as Interacting Particle Systems^^: 

One might think of the parking processes P„,n > described in definition 12.11 
as discrete-time versions (i.e. skeletons) of "finite" interacting particle systems 
{'qf'^)t>o with initial condition ryo = 0^ and flip rates cix, rf) given by 

c(x,77) = • l{e,(^)|A('')e5<^'} Vxe Z^Tye {0,1}^' 

It is not difficult to see that this is indeed the case, provided 5^,'^'' is a decreasing set 
(The foregoing correspondence being not correct otherwise, since non-decreasing 
parking schemes generate non-markovian parking processes.) 

In this setting it would be natural to define the properly infinite interacting 
particle system {r]t)t>Q, defined by the flip rates 

cix, ri) = l{e,(^)|A(-)es(^)} Va; e 77 G {0, 1}^' 

and same initial condition r/o — 0^ and construct (i.e. couple) {r|["'^)t>{)^n > and 
('7i)t>o together making use of the same Poisson marks to conclude that 

lim ?7f'-' r/^^ r/oo hm r/t a.s. 

t — >oo t — >oo 

Since the sequence {ri^^)n>o is probabilisticaly indistinguishable from the the se- 
quence of jamming limits {Pn)n>o from the foregoing sections (i.e. {r]^x^)n>a ^ 
(_P„)„>o), we may interpret rjac as P, "the limit jamming limit" of section |21 (defi- 
nition 

All this reasoning is meaningful and in fact it proves in straghtfoward fashion the 
existence of the limit measure /i (at least for the markovian case). Yet one should 
bear in mind that the proofs carried out in sections 0] and [S] rely on the specifics of 
section|31 Put it in other terms, the sole existence of p is not enough to accomplish 
the results of this paper. 



^"'^It is shown therein that E(|P„|) ~ 0,4324 ■ 2n for the one-dimensional nearst neighbour 
exclusion scheme we touched on above, where \Pn\ is short for g^(ti) Pn{x) as in section [Tl In 
fact, Mackenzie thought in terms of unit intervals instead of sites/particles and made calculations 
for the expected vacant space in the jamming limit instead of the expected number of occupied 
sites. Here we have translated his terminology to our current language. 

■^^We refer the interested reader to R. Durrett's review article 1' for a quick, yet meaningful, 
introduction to this subject. 

^■^As we already paid attention to, we are abusing terminology and using the jamming limit 
Pn to denote the underlying parking process. 

■^"^i.e. for ^1,^2 S S'f,'^', if ^1 =4 £.2, then ^2 S sif ^ => 5i G si'^^; where '=<;' is taken for the usual 
partial order of spin configurations. 
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Parking Processes with Boundary Conditions Other than Null: 

It is easy to see that whatever boundary condition we use in definition l2 . II we always 
obtain the same thermodynamic jamming limit P of definition 13. 61 More precisely 
if the sequence (-Pn)n>i is defined in the same probability space of the former 
sequence {Pn)n>i'i exactly as prescribed in definition 12.11 with the sole exception 
that we modify step(i) to 

Step(i) F„|AW ^ 0"^''!" ; P„ \ (z'^ \ Ai^^) ^ (z"^ \ A^f)) 

for an arbitrary rj G {0, 1}^ , then 

Pn ^ P ^ P71 cl.S. 

To see this we only have to modify lemma slightly as follows: 

For each ui £ {0, 1}^ and x G Z'^, i/ A„ D A {{x}) (uj), then Pn+u and Pa({x}) ci^s 
identical over A {{x}); that is Pn+u\A{{x}) = Pa{{x})\ -^{{x})', 

and to take n sufficiently large in theorem 13.31 so that A{Am^) C A^fL^- 
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